DISCRETE DIRAC OPERATORS, CRITICAL EMBEDDINGS AND IHARA-SELBERG 

FUNCTIONS 



MARTIN LOEBL AND PETR SOMBERG 

Abstract. The aim of this paper is to formulate a discrete analog of the claim made by Alvarez-Gaume et 
al., ([a]), stating that the partition function of the free fermion on a closed Ricmann surface of genus g is 
certain linear combination of 2 2fl Pfaffians of Dirac operators. Let G = (V, E) be a finite graph embedded in 
a closed Riemann surface X of genus g, x e the collection of independent variables associated with each edge 
e of G (collected in one vector variable x) and £ the set of all 2 2ff Spin-structures on X. We introduce 2 2 9 
rotations rot B and (2\E\ X 2\E\) matrices A(s)(x), s £ E, of the transitions between the oriented edges of G 
determined by rotations rot s . We show that the generating function of the even sets of edges of G, i.e., the 
Ising partition function, is a linear combination of the square roots of 2 2ff Ihara-Selbcrg functions I (A(s)(x)) 
also called Feynman functions. By a result of Foata and Zeilberger I(A(s)(x)) = det(7 — A' (s)(x)), where 
A'(s)(x) is obtained from A(s)(x) by replacing some entries by 0. Each Feynman function is thus computable 
in a polynomial time. We suggest that in the case of critical embedding and bipartite graph G, the Feynman 
functions provide suitable discrete analogues of the Pfaffians of Dirac operators. 



1. Introduction 

It is well known, see e.g. [mj, [k] and references therein, how to formulate critical embedding of a finite 
graph in a closed Riemann surface in such a way that one can read off from the collection of angles of 
the embedding, termed discrete conformal structure, the critical values of independent variables (i.e. the 
edge weights or coupling constants in the physical terminology) attached to the edges of the Dimer and the 
Ising problems on G. It is rather attractive task to study whether some of the properties associated to the 
notion of criticality in statistical physics may also be derived from the collection of data attached to critical 
embeddings. 

The main theme of this paper is the formulation of certain discrete analog of the claim made by Alvarez- 
Gaume et all, see [aj, that the partition function of free fermion on a closed Riemann surface of genus g is 
a linear combination of 2 2g Pfaffians of Dirac operators. The theory of free fermion is generally accepted to 
be closely related to the criticality of both the Dimer and the Ising problems on graphs. 

Quite recently the Dimer problem and its determinant type solution received considerable interest, see e.g. 
[k], [cj. The authors of [bt] study the critical Ising model by reducing it via the determinant type method 
to the Dimer model. The case of planar graphs is well-understood in this setting, [kj. However, as proved in 
[c] , if one wants to obtain a discrete analog of the claim of |a] for the Dimer model in a surface of a positive 
genus, one has to take into account global restrictions on the graph embedded in the surface. 

For the sake of completeness, we review in the Appendix the combinatorial approach used to describe de- 
terminant type reduction of the critical Dimer model on Riemann surfaces of positive genus. The obstruction 
emerging in this description is centered around the notion of Kasteleyn flatness. 

In this paper we propose to overcome the previously mentioned limitations and obstacles for the determinant- 
type reasoning by replacing the determinants by so called Ihara-Selberg functions of the graph G. This in 
turn allows to rewrite the generating function of the even sets of edges of G (the Ising model partition func- 
tion associated to the graph G) as a linear combination of 2 2g square roots of Ihara-Selberg functions which 
we call Feynman functions. We build the Feynman functions in a combinatorial way in order to capture 
the emergence of the rotations rot s for Spin structures s in the analysis of the Ising partition function, and 
realization of rot s via the quadratic forms on Hi(X, F2). This approach makes use of the original treatment 
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of the Ising partition function via Ihara-Selberg functions by Sherman, see [s], pj. Our basic notion of a 
g— graph G 9 allows convenient treatment of a planar model for the embedding of the graph G in a closed 
Riemann surface of genus g. The g— graphs have been successfully used recently in [lmj in a related situation. 
Moreover, the planar model may improve the insight into the analysis of the limiting process, which has not 
been done yet. The rotations are treated in a different (and perhaps more elegant) way by Johnson Rl. 

Finally, we show that the discrete Dirac operators appear very naturally in the context of the Feynman 
functions and comment on the relationship with continuous counterpart called analytic torsion, see |rsl) . |rs2] , 

1.1. Critical embeddings of graphs into closed Riemann surfaces. For the material reviewed in this 
subsection, see e.g. [mj and references therein. Let us consider a triple (X, G, tp), where X is a closed 
Riemann surface, G — (V, E) a graph and ip : G X an embedding. The image of ip defines a CW 
decomposition of X and X \ ip(G) is a disjoint union of open faces. 

The metric realization of X, compatible with an atlas {Ui, (pi : Ui — »■ C| |Jj Ui = X} covering X, looks as 
follows. For a finite set of points {Pj}j on X there exists a flat metric with conical singularities at {Pj}j 

such that the cone angles 6p j fulfill the Gauss-Bonnet formula x{X) — YljO- — air")' an< ^ there is no other 
restriction on them. Because the Euler characteristic xPO is even, one can choose 9p j to be odd multiples of 
2ir for each j, for example. Recall that the local model for a neighborhood of a conical singularity in Pj is the 
standard cone C{0) := {(r,j3)\r £ G R/0R}/{(O,/3) ~ (0,/3')} with the metric g C (e) ■= (dr) 2 +r 2 (d/3) 2 . 

An embedding of a graph G in a surface X induces the dual graph G* of the embedding. We regard 
G* as an abstract graph with natural embedding into X, so that each vertex of G* lies on the face of the 
embedding of G it represents. The central notion related to the couple G, G* is of diamond graph - for 
simultaneous embedding of G and G*, the diamond graph G + has the vertex set equal to V(G) U V(G*) and 
the edges connecting the end- vertices of each dual pair of edges e, e* into a facial cycle F(e) of G + , which is 
a 4— gon called a diamond (see Figure [TJ. 
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Figure 1. Diamond F(e) of edge e. 

We say that an embedding of G is critical if it is isoradial, i.e., in any local trivialization realized by 
complex plane is each face of G a cyclic polygon (the face is possible to inscribe into a circle.) Moreover, the 
circumcenters of the facial polygons are contained in the closures of faces and all the facial polygons have 
the same radius. An equivalent description based on diamond graph G + looks as follows. The embedding 
cp of G + is critical if each of its faces F(e), e G E, is a rhombus, i.e., the following conditions hold true (for 
all i) with respect to the induced conformal class of metrics on <p(X): 

(1) The diagonals of each rhombus (in Im(cpi o ip)) are perpendicular, 

(2) The lengths of sides of all rhombi (in Im(ipi o <p)) are the same. 

In particular, the first property is independent of the choice of trivialization, i.e. of the index i together 
with Im(tpi o ip), because the transition maps ipi o ipj 1 are conformal and so angle preserving. The second 
condition already depends on the choice of a metric in a given conformal class of metrics. 
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1.2. Dimer and Ising partition functions. We associate an independent variable x e with each edge 
e € E of G. A subset E' C E ai edges is called perfect matching or dimer arrangement, if the induced graph 
(V, E') has each vertex of degree one. Let V(G) denote the set of the perfect matchings of G. We define the 
dimer partition function of G by 

(1) v(g,x)= y, n^> 

MeV(G) eGM 

where x = {x e ) e ^E is the vector of edge weights. 

A subset E' C E of edges is called even if the induced graph (V, E') has each vertex of even degree. We 
denote by £{G) the set of even subsets of edges of G. 

The generating function of the even sets of edges of G is defined by 

(2) £{G,x) = I]>- 



It is well known (see e.g. |lm| ) that £{G,x) is equivalent to the Ising partition function on G defined by 



(3) 4 SinS (/3) == 4 SinS (*) 



°f 3 ° VeeE 



where J e (e G E) are the weights (coupling constants) associated with edges of the graph G and (3 the scale 
(inverse temperature), and 



(4) z™*{ X )= £ n ' 



o(u)cr{v) 

cr:V->{l,-l} e={tt,«}£B 



1.3. The aim of the article and its main results. We denote by E the set of Spin structures on Riemann 
surface X . Following the articles [!],[]], |lmj . we associate the rotation to each spin structure. This enables 
us to define (2\E\ x 2\E\) matrices A(s)(x), s £ E in a way that each A(s)(x) is the matrix of the transitions 
between the oriented edges determined by the rotation rot s corresponding to the spin structure s. The main 
result of the article is 

Theorem 1. Let G be a graph embedded in a closed Riemann surface X of genus g and let E be the set of 
Spin structures on X. Let us denote for each s £ E by A(s)(x) the matrix introduced in Definition \2.10\ 
Let F(A(s)(x)) be the Feynman function defined in Definition H77H fl s the square root of the Ihara-Selberg 
function I(A(s)(x)). Then £(G,x) is a linear combination of 2 2g functions F(A(s)(x)), s £ E. Moreover, 
each F(A(s)(x)) is computable in a polynomial time. 

Furthermore, let G be a bipartite graph critically embedded in the Riemann surface X of genus g equipped 
with a flat metric with conical singularities in such a way that no conical singularity is located at a vertex of 
the graph G. Then for each Spin structure s € E the matrix A2(s) introduced in DeHnition \2. lb\ is related 
to the discrete Dirac operator associated to s (see [c] for its definition) in a way described in Corollary 1 2. 19\ 



Acknowledgements. The authors would like to thank David Cimasoni for extensive discussions about 
discrete Dirac operators and Gregor Masbaum for extensive discussions and suggestions about rotations and 
spin structures. 



2. Ihara-Selberg functions and Feynman functions 

In this section we assume the graph G is embedded in a closed Riemann surface X of genus g. We suggest 
to consider the Ising partition function £{G,x) instead of the Dimer partition function V(G,x), and the 
square root of certain Ihara-Selberg functions on the graph G (which we call the Feynman functions) instead 
of the determinant. Similar but much less advanced structure can be found e.g. in jsrn] under the notion of 
Ising preholomorphic observable. 

3 



2.1. Ihara-Selberg function. Let G = (V,E) be a graph. For e£fi we denote by o e an orientation of e, 
and o~ x the reversed directed edge to o e . As above, let x — (x e ) e ^E be the formal variables associated with 
edges of G. If o is any orientation of the edge e, we associate the new variable x Q with it and always let 

X Xe~ 

We consider an equivalence relation on the set of finite-length sequences (zi, . . . , z n ) satisfying z\ — z n : 
each such sequence is equivalent with each of its cyclic shifts. The equivalence classes will be called circular 
sequences. 

A circular sequence p — Vi, Oi, V2, 02, o n , v n +i with v n+ i — v\ is called a prime reduced cycle if the 
following conditions are satisfied: Oi <E \p e ,o~ x : e S E}, Oi oL^ and (01, ...,o n ) ^ Z m for some sequence 
Z and a natural number m > 1. We say that the ordered pair (oi,Oi + \) is a transition of p at i>j+i, and 
(o n , Oi) is a transition of p at t>i. We denote by p _1 the prime reduced cycle which is the inverse of p. 

Definition 2.1. Let G = (V,E) be a graph and assume that the vertex set E is linearly ordered. Let M 
be 2\E\ x 2\E\ matrix with entries m(o,o'), 0,0' s {oe^o^ 1 : e 6 where we think of m(o, o') as the 
weight of transition between directed edges o, d of G. If p is a prime reduced cycle then we let M(p) = 

Yi(o,o') a transition of p TO (°7 ° )• 

We denote the set of prime reduced cycles of G by Q. The Ihara-Selberg function associated to G is 

(5) /(M)=IJ (1-^(7)) 

where the infinite product is defined by the formal power series: 

(6) II( 1 - M (7))=E(- 1 )^ II M (7) 

and the sum is over all finite subsets J- of Q. 

A theorem of Foata, Zeilberger (see [fe]), generalizing the seminal theorem of Bass (see [b]), states: 
Theorem 2. 

(7) I(M) = det(I - M'), 

where M 1 is the matrix obtained from M by letting m'(o, o') = if d = o _1 and m'(o,o') = m(o, d) 
otherwise. 

2.2. Combinatorial model for closed Riemann surfaces of genus g. In this section we restrict our- 
selves to the following standard representation of closed Riemann surface X of genus g: we regard A as a 
regular Ag— gon R (called the base polygon) in the plane with sides denoted anti-clockwise by z\, . . . , Z4 g , 
and the pairs of sides Zi, z~^ 2 and 2j+i, z^r 3 , i — 1, 5, . . . , i(g — 1) + 1, identified. This defines a flat metric 
on X with one conical singularity of angle 2n(2g — 1). We only consider embeddings of graphs on X which 
meet the boundary of R transversally. Moreover, to simplify the arguments, we only consider embeddings 
of graphs where each edge is represented by a straight line on X . The general embeddings may be treated 
analogously. 

We now describe how an embedding of a graph in a Riemann surface can be used to make its planar 
drawing of a special kind. We follow [I], [lmj . 

Definition 2.2. The highway surface S g consists of the base polygon R and the bridges R\, . . . , i?2g, where 

• Each odd bridge i?2i-i is a rectangle with vertices x(i, 1), . . . ,x(i,4) numbered anti-clockwise. The 
bridge is glued to R so that its edge [x(i, l),x(i, 2)] is identified with the edge [^4(i-i)+i, ^4(i-i)+2] 
and the edge [x(i,3),x(i,A)] is identified with the edge [z4(i-i)+3, Z4U-i)+4\- 

• Each even bridge i?2i is a rectangle with vertices y{i, 1), . . . , y{i, 4) numbered anti-clockwise. It is 
glued with R so that its edge [y{i, l),y(i,2)] is identified with the edge [24(^-1)4-2, ^4(i-i)+3] and the 
edge [y(i, 3), y{i, 4)] is identified with the edge [2:4(^1)4-4, 24(^1)45] (the indexes are always considered 
modulo Ag.) 

There is an orientation-preserving immersion $ of S g into the plane which is injective except that for each 
i = 1, . . .g, the images of the bridges i?2i and i?2i~i intersect in a square (see Figure [5]). 
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Now assume the graph G is straight-lines-embedded into a closed Riemann surface X of genus g. We 
realize X as the union of S g and an additional disk i?oo glued to the boundary of S g . By an isotopy of the 
embedding we may assume that G does not meet the disk i?oo and moreover, all vertices of G lie in the 
interior of R. We may also assume that the intersection of G with any of the rectangular bridges Ri consists 
of disjoint straight lines connecting the two sides of Ri glued to the base polygon R. The composition of the 
embedding of G into S g with the immersion $ yields a drawing tp of G in the plane, where each edge of G 
is represented by a piece- wise linear curve (see Figure [2]). A planar drawing of G obtained in this way will 
be called g— graph and denoted by G 9 . Observe that double points of a g— graph can only come from the 
intersection of the images of bridges under the immersion <1> of S g into the plane. Thus every double point 
of a special drawing lies in one of the squares $(i?2i) H $(i?2i-i)- 

Definition 2.3. Let G be embedded in S g and let e be an edge of G. By definition, the embedding of e 
intersects each bridge Ri in disjoint straight lines. The number of these lines is denoted by 7-j(e). For a set 
A of edges of G we let r(A) be the vector of length 2g defined by r(A)i = J2 e eA r i( e )- 

Definition 2.4. If p is a prime reduced cycle of G then we denote by p 9 the image of p in G 9 . 

Clearly, p — s> p 9 gives bijective correspondence between the prime reduced cycles of G and the prime 
reduced cycles of G 9 . 



Figure 2. Immersions of edges {a, b}, {c, d} and {e, /} which cross a side of R. 

2.3. Rotations. Let G be a graph embedded in a closed Riemann surface X of genus g and let G 9 be its 
g— graph. We recall that each edge of G in G 9 is represented by a piecewise linear curve. Let p be a prime 
reduced cycle of G and let p 9 the corresponding prime reduced cycle of G 9 . We will introduce 4 s rotations 
ofp 3 . 

First of all, we denote by the 0— vector of length 2g and in analogy with the usual definition of the 
rotation of a regular closed curve in the plane we define roto by setting roto(p 9 ) — J2t Vo (*) (mod 2), where 
we sum over all the transitions of the linear parts of p 9 . If the transition t consists in passing from directed 
segment e to directed segment e' then yo(t) = Zo(t)(2Tr)~ 1 , where zo(t) is the angle of the transition. The 
angle zo(t) is negative if the transition is clockwise, and zn(i) is positive if the transition is anti-clockwise 
(see Figure [3|). Consequently, for each arithmetic vector s G F% 9 (F% is the field with two elements) we let 
y s {t) = yo(t) + sr({e}) and rot s (p 9 ) — ^ t y s {t) — sr(p) + rot (p 9 ) (mod 2). Here sr denotes the scalar 
product of vectors s, r. Observe that for each s, (— l)' ot »(p 9 ) — (—iy ot s((p s ) x ) _ 

Definition 2.5. We introduce the equivalence relation on the set of prime reduced cycles of G: we say that 
Pi is equivalent to P2 if pi = {P2)~ l ■ The set of equivalence classes of this equivalence relation is denoted by 
[Q] . Analogously, we introduce the equivalence relation on the set of prime reduced cycles of G 9 : we say that 




e 
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pf is equivalent to p 9 2 if pf = • The set of equivalence classes of this equivalence relation is denoted by 

[G] g . 



V V 
+tt/2 -tt/2 

Figure 3. Angle of the transition. 

The previous considerations allow to introduce the Feynman functions: 

(8) F( G ,x,s)=j2(-^ m n (-ir M7) n-- 

[F\a [ 7 ]e[f]» ee 7 

where the sum is over all finite subsets [J 7 ] 9 of [Q] 9 . The following theorem appears in [I]. We assign to 
s e F% 9 the sign according to sign(s) = (-l)rii=i,3,...2 9 -i s 4 s j+ ^ 

Theorem 3. Let G be a graph with each vertex of degree equal to 2 or 4, embedded into a closed Riemann 
surface X of genus g. Then 

(9) £(G,x) =2- 9 ^2sign(s)F(G,x,s). 

Theorem [3] can easily be extended to general graphs. 
Theorem 4. Let G be a graph embedded into a closed Riemann surface X of genus g. Then 

(10) £(G,x) = 2~ 9 siga(s)F(G,x,s). 

Proof. We will construct a graph G' along with its embedding into X so that the degree of each vertex of G' 
is equal to 2 or 4, such that there are two subsets Z, O of E{G') and a bijection / : E{G') \ (ZUO) — > E{G) 
inducing 

£(G, x) = £(G' , z)\ Ze ;— Xf({ ^ if e(ZUO;z a :=0 if e<EZ;z c : = l if eEO 



and, for each s G F? 9 



/(=) 

2 ' 

F(G, a;, s) = F(G , Z, s)| Zs ~ x ^ (e) if e ^ZuO;z c :=0 if eeZ;z c : = l if eeO- 

Theorem |4] then follows from Theorem[3] We construct graph G' in two steps. Let OD denote the set of the 
vertices of G of an odd degree. We start with O = Z = 0. 

Step 1. If OD then it is a standard observation of the graph theory that G has a set of edge-disjoint 
paths so that each vertex of OD is an end-vertex of exactly one of the paths, and all the end- vertices of the 
paths are among the elements of OD . In particular, \OD\ is even. Let P denote the set of the edges of 
these paths. We construct graph G\ from G by adding, for each edge {u, v} £ P, a path of length 3 with 
end- vertices u, v (see Figure 0]). We further let Z be the set of all the edges of G\ \ G. Note that G\ has all 
degrees even. 

Step 2. If v is a vertex of G\ of an even degree bigger than 4 then we modify G\ by splitting the degree 
of v by introducing the new splitting edge; the operation can be read off the Figure 03 We repeat this step 
until the resulting graph G' = (V, E') has all degrees equal to 2 or 4. Finally we let / be the tautological 
injection of E into E' . 

It is straightforward to realize that all assumptions of the construction are satisfied after application of 
finite number of these steps. 

□ 
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Figure 5. 0:=0U{ the new splitting edge } 

In order to associate rotations to quadratic forms, we first need to study self-intersections of the prime 
reduced cycles which traverse each edge of the graph at most once. 

2.4. Self-intersections of prime reduced cycles. Let p be a prime reduced cycle of G which traverses 
each edge of G at most once. For each vertex v of G, let p(v) denote the set of the directed edges of p 
incident with v, and let P(p, v) denote the partition of p(v) into pairs which correspond to the transitions 
of p at v. If a prime reduced cycle p traverses each edge of G at most once then the transitions of p at v are 
well described by the directed chord diagram diag(p,v) (see Figure[6]): 

Definition 2.6. Let p be a prime reduced cycle of G which traverses each edge of G at most once. The 
directed chord diagram diag(p, v) is obtained by taking the cyclic ordering of the edges of G incident with 
v and induced from the embedding of G in X, and by introducing the directed chord (e, e') for each class of 
P(p, v) consisting of an orientation of e followed by an orientation of e'. 

We define the number of self-intersections of p as the number of the pairs of intersecting chords of 
diag(p, v), v G V. 




d 




a 



Figure 6. Directed chord diagram diag(p, v); p — ac . . . db . . . a 



2.5. Quadratic forms. A subset of edges E' C E of a graph G embedded in a closed Riemann surface of 
genus g is called even if each degree of the graph (V, E') is even. The set of even subsets of G is denoted by 
£(G). Let H := Hi(X, F2) be the first homology group of X with coefficients in the field F2. We observe that 
two even subsets of edges A, B belong to the same homology class in H if and only if r(A) = r{B) (mod 2) 
(see Definition 12.31 for the definition of r). We have a basis a\, b\, . . . , a g , b g of H, where a*, &j correspond to 
the even subsets Ai,Bi satisfying r(^4»)2i-i = 1 and r(Aj) = otherwise, and r(Bi)2i — 1 and r{Bj) = 
otherwise. 

Recall that H carries a non-degenerate skew-symmetric bilinear form called (mod 2) intersection form, 
denoted by In the basis chosen above, it is given by a, -a,j —bi-bj = 0andaj-6j = 8ij for alH,j = 1, . . . , g. 

Definition 2.7. A quadratic form on (H, ■) associated to the bilinear form '•' is a function q : H F2 
fulfilling q(x + y) = q(x) + q(y) + x ■ y for all (x, y G H). 

We denote the set of quadratic forms on H over F2 by Q. The cardinality of Q is 4 9 . 

Let E' G £(G) be an even subset of edges of G. For each vertex u, let E'(v) denote the subset of edges 
of E' incident with v, and given the embedding of G in X, let N(E',v) denote a partition of E'{v) into 
pairs which correspond to the non-crossing undirected transitions at v on X (see Figure EJ. Notice that 
N(E', v) need not be unique. We say that a prime reduced cycle p of G is consistent with U t , 6 vA^(£", v) if 
each transition of p is a directed transition of an element of U^ e yiV(i? / , v). 



Figure 7. If £'(«) = {a, b, c, d, e, /} then JV(£' , v) = af, be, cd 

Let us recall Definition 12.51 of the equivalence on the set of prime reduced cycles of G. 

Observation 2.8. The set of all transitions {J vl zy N [E' , v) uniquely determines the set W(E') = {[pi], ■ ■ ■ , [pk]} 
of equivalence classes of edge-disjoint prime reduced cycles which are consistent with U V £vN(E' , v), and 
UiPi = E' . Moreover, for each i we have rotoipf) = 1+ ^ e number of the self-intersections of pf (see 
Definition \2.6\) . 

Proof. The first part is straightforward. The second part is clear since no edge is traversed more than once 
by any of the Pi's. 



We recall that £ (G) denotes the set of the even subsets of edges of graph G. Let H(G) C H denote the 
subgroup consisting of the homology classes realized by elements of £{G). 

Definition 2.9. For each s G F^ 9 , we define a function q' s : £{G) — > F2 by 




□ 




pew(E') 



where p a is the realization of p in the g— graph G 9 (see Definition 12 .4j) . 
Theorem 5. For each s G F^ 9 , the function q' s induces a function 
(12) q s : H(G) -> F 2 
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by 

(13) q s ([E'])= q ' s (E'). 
Moreover, each q s is a quadratic form restricted to H(G). 

Proof. We first show that qo is a quadratic form. Let W(E') — {[pi], . . . , [pk]} be as in Observation ^. 81 We 
have, for each i, roto(pf) = 1+ the number of the self-intersections of pf (Observation 12. 8[ ). Moreover, by 
the choice of \J V £vN(E', v), the self-intersections of each pf as well as the crossings of pf with pP, i ^ j, 
lie outside of the base polygon R. We also recall that two closed curves in the plane have always an even 
number of crossings. 

Hence, q'o(E') — J2i=i r (^')'2i-i r {E')2i (mod 2). This means that, by Definition 12.71 go coincides on 
H(G) C H with the quadratic form on H whose value on each of the basis vectors a, and bi is zero. 

Furthermore, rot s (p 9 ) = sr(p) + roto(p 9 ) (mod 2) and thus q' s (E') = sr(E') + q' {E') (mod 2). Hence 
each q s is well-defined and coincides on H{G) C H with a quadratic form. 

□ 

We recall that there is a natural bijection between the set £ of the Spin structures of X and the quadratic 
forms on the F%— valued first homology classes of X (see []]). Hence Theorem [S] associates a Spin structure 
to each rot s , s 6 F 2 9 , and this map is bijective. From now on we will consider the rotations rot s indexed by 
Spin structures on X. However, we still have rot s determined by the derived drawing of G realized by the 
g— graph G 9 . 

The usefulness of quadratic forms in studying rotations rot s was suggested to us by G. Masbaum, who 
also suggested Theorem [6] below. 

Theorem 6. Let p be a prime reduced cycle of G so that each edge of G is traversed at most once in p, and 
let s € S. Then (mod 2) 

(14) rot s {p 9 ) = 1 + the number of the self-intersections of p + q s (p), 
where p 9 is the realization of p in the g— graph G 9 . 

Proof. By the definition of rot s it suffices to prove the statement for s = 0. We have (mod 2) 

roto(p 9 ) = 1 + the number of the self-intersections of p 9 = 

a 

1 + the number of the self-intersections of p + f(p)2i-if(p)2i = 

t=i 

1 + the number of the self-intersections of p + <jfo(p)- 

□ 

2.6. Feynman functions. 

Definition 2.10. Let G be a graph embedded in a closed Riemann surface X of genus g. Let E° = U{o e , o^ 1 : 
e e E}, hence \E°\ = 2\E\. To each Spin structure s £ S we associate A(s)(x), the \E°\ x |_E°|-matrix with 
entries d s (o,o') = (-l)»»(°' ')(-l)"( )a£, where 

(1) k(o) = if o is contained in the interior of R, 

(2) k(o) = —3/4 if the segment of the embedding of o outside R in G 9 is oriented oppositely to the 
anti-clockwise orientation of the boundary of R (as directed edge (a, b) in Figure [5]), 

(3) k(o) — 3/4 if the segment of the embedding of o outside R in G 9 is oriented in agreement with the 
ant i- clockwise orientation of the boundary of R (as directed edge (a, b) in Figure [2]). 

We further define A' (s)(x) by declaring d' s (o, o') =0ifo' = o _1 and d' s (o,o') — d s (o, o') otherwise in A(s)(i). 

We are going to introduce the Feynman functions, see Definition 12.121 below. As its first ingredient we 
need to extend Definition 12.61 to the number of self-intersections of general prime reduced cycles, i.e., the 
prime reduced cycles that can go through an edge more than once. This can be done for instance in the 
following way. We define the infinite graph G by replacing each edge e by an infinite sequence of edges 
ei, • • • , ej, • • ■ with the same end- vertex as e. We embed G in X so that we 'thicken' the embedding of each 
edge e of G, and embed the edges ei, • • • , e^, • • • to this thickened part of e so that they are piece- wise linear 
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and internally disjoint. Next, for each prime reduced cycle p of G whose circular sequence of directed edges 
is (ai . . . , dfc), a,j being an orientation of edge e(j) of G (possibly e(j) — e(l) for j ^ I), we define prime 
reduced cycle p in G by replacing each aj by the same orientation of e(j)j- It is important that the prime 
reduced cycle p uses each edge of G at most once. We thus define the number of the self-intersections of a 
prime reduced cycle p of G as the number of the self-intersections of the prime reduced cycle p of G (see 
Definition ^. 6[) . The number of the self-intersections of a prime reduced cycle p 9 of G 9 is defined analogously. 

We note that for the prime reduced cycles of G containing each edge of G at most once this is consistent 
with Definition 12.61 We also note that the following basic property is satisfied. 

Observation 2.11. Let p be a prime reduced cycle of G. Then we have (mod 2) 

1 + the number of the self-intersections of p 9 = roto(p 9 )- 

Definition 2.12. Let s € £ be a Spin structure. We define 

(15) F(A(s)(x)) := V J[ (_1)9»(7)+ the number of self-intersections of 7 TT ^ 

[F\ [ 7 ]e[.F] ee 7 

Proof, (of the first part of Theorem [TJ 

We note that 

F(A.(s)(x)) = S V — l)!^! J^J ^_^^l+9s(7)+ the number of the self-intersections of 7 II x _ 

[7] [7]e[^] e G7 

= ^(-1)1^1 J] (-1)^11^ 

Hence F(A(s)(x)) is the square root of I{A(s)(x)). By Theorem 7(A(s)(x)) = det(7 - A'(s)(x)) 
and each F(A(s)(x)) is computable in polynomial time. The first part of Theorem [1] is now a corollary of 
Theorem 0J 

□ 

2.7. The Arf-invariant formula. The linear combination appearing in Theorem [T] may be understood as 
the Arf-invariant formula 

(16) £(G, x) = 2~" ^(-l) Arf ^)F(A( S )( a ;)). 

This formula was suggested by G. Masbaum and independently by D. Cimasoni, [cc\ A proof following [gj, 
[I] is work in progress (see [lmmj ) . We only indicate here the emergence of £ (G, x) in the right- hand- side of 
the formula above. We know by Theorem [5] that for each even set E' € £{G) and for each seS, 

(17) q s (E')= ]T (l + rot s (p 9 )) (mod 2). 

Next, the basic result about the Arf invariant (see e.g. Lemma 2.10 of |lmj ) states that 

(18) 1 = 2- a J2(- 1 ) Ali{9s) (- 1 ) 9s(E ' ) - 

sGS 

Hence, for each E' 6 £(G), YieeE' Xe a PP ears m the expansion of the right-hand-side of the Arf-invariant 
formula. It is needed to show that the remaining terms in the expansion are equal to zero. 

2.8. Discrete Dirac operators. So far we considered an arbitrary embedding of a graph G in a Riemann 
surface X. Now, let G — (W,B,E) be bipartite and the embedding critical. The vertices of W are called 
white and the vertices of B are called black. We also assume in this subsection that the conical singularities 
are not located at the vertices of G. 

We recall the notation E° = Ll-fo^o^ 1 : e £ E}. Let W° be the subset of E° consisting of the edges 
directed from its black vertices to its white vertices. Analogously, let B° — E° \ W° be the set of edges 
directed from its white vertex to its black vertex. 

The key construction has the following structure: 
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Definition 2.13. Let Tq be the directed transition graph of the orientations of edges of G, i.e. V(Tg) = E° 
and (o, o') € E(Tq) if the head of o is the tail of d . 

Observation 2.14. The graph Tq is a directed bipartite graph, Tq — (W°,B°,E(Tq)) ) and the matrix 
A(s)(x) is a weighted adjacency matrix of To. 

The next observation follows directly from Definition 12.101 

Observation 2.15. Let w° G W° be an orientation of the edge e entering the vertex w 6 W . Let b° be a di- 
rected edge leaving w and entering vertex b G B. Then the entry A(s)(x)(w° , b°) equals 7(s, w°)e la ^ w ,b ^Xb a! 
where a(w°, b°) = l/2zo(w°, b°) is half of the angle of the transition from w° to b°; j(w ) equals a complex 
number depending only on s and w° . 

Let l(e*) denote the length of the dual edge e* of an edge e of G. 

Definition 2.16. We denote by Aa(s) the matrix obtained by taking the square of each entry of A(s)(x) and 
substituting the vector of the lengths of the dual edges for the squares of the variables: (x e ) 2 := l(e*), e G E° . 

Corollary 2.17. (of Observation ^. 15\) IfuP^w 1 are two directed edges of E° entering the same vertex then 
the row of A 2 (s) indexed by w° is a complex multiple of the row o/A 2 (s) indexed by w . 

The discrete Dirac operator T)(s) of a bipartite graph G = (V, E), V = W L) B, corresponding to a Spin 
structure s G E, is defined as the weighted adjacency (\V\ x \V\)— matrix; [k. contains the definition for 
the planar graphs and Jc] defines half of (or, the chiral part of) T>(s), which uniquely determines T>(s). We 
observe in this subsection that the matrix A 2 (s) is closely related to T>(s). 

To that aim we fix the coordinate chart on the Riemann surface such that a black vertex v G B lies 
at £ C in this local chart. The directed edges emanating from a white vertex and terminating at v are 
denoted e\, . . . , e%.. Similarly, we fix an edge e v w emanating from v and terminating at some white vertex 
w G W. The (i!w;)-matrix element of the Dirac operator is Di vw \ = Z((e^,)*)e !Q , where a is the angle between 
the edge e v w and the real axis of the local coordinate chart. Note that the change in local coordinate chart 
does not preserve the real axis, but does preserve the angle a. On the other hand the (e", e^)-entry of 
our matrix A 2 (s) is equal to Z((e^)*)e lQi , where ai is the angle between e^,e^ measured anticlockwise, see 
Figure [3] In particular, in the local coordinate chart in which the edge e\ lies on the real axis, we have 
a = ai. Equivalently, the matrix coefficients in a given row of our matrix A2(s) indexed by directed edge e 
are (constant) multiples of the row in the Dirac matrix corresponding to the terminal vertex of e. The same 
applies when one starts with the white vertex instead of the black one. 

The basic observation about the matrices A 2 (s) and T)(s) can be formulated in the following way. 

Corollary 2.18. Let w° G W° be an orientation of the edge e entering the vertex w G W. Let b° be a 

directed edge leaving w and entering vetex b G B. Then A2(s)(w° ,b°) = c(w°)T>(s)(w,b) , where c{w°) is a 
complex number depending only on w° . 

The last part of Theorem Q] is contained in the next Corollary. 

Corollary 2.19. Each matrix A 2 (s), s G S a spin structure, may be obtained from the discrete Dirac 
operator T>(s), whose chiral part is defined in [gj, by finite number of operations: 

(1) Adding identical copy of a row, 

(2) Multiplying a row by a complex number, 

(3) Adding 2\E\ — \V\ zero-entries to each row. 

Remark 2.20. The properties that graph G is bipartite and critically embedded, and the condition that the 
conic singularities of the metric are not located at the vertices of G, are needed in [nj to show that constant 
functions belong to the kernel of the discrete Dirac operators. 

2.9. Fermionic expansion of the Feynman functions. In the last part of the paper we rewrite each 
Feynman function as the alternating sum of the traces of the skew-symmetric powers of matrices A'(s)(x). 
The skew-symmetric powers are closely related to the analytic torsion on Riemann surfaces. These consid- 
erations may lead to a realization of the analytic torsion as a continuous limit of the Feynman functions. 

Let T G be the directed graph obtained from Tq by deleting all directed edges of the type (o, o" 1 ), where 
o is an orientation of an edge of G. Clearly, A'(s)(x) is the weighted adjacency matrix of the directed graph 
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T' G . We recall that each Feynman function F(A(s)(x)) is the square root dct(7 — A'(s)(x)) (Theorem [2]) . 
We first describe the directed graph T' G in terms of the graph G. 

The set of the vertices of T' G is W° U B° . Let e = {b, w} be an edge of G between black vertex b and white 
vertex w. The edge e gives rise to two vertices of T' G , namely — (w,b) and e w = (b,w). If the graph G 
is critically embedded in a closed Riemann surface X , the vertex e;, of T' G is represented by a point on the 
embedding of e 'near' the vertex 6, and the vertex e w of T' G is represented by a point on the embedding of e 
'near' the vertex w. (see Figure 

Let toi, . . . , Wk be the white vertices incident to b in G and different from uu, and let e l = {b, Wi}. Then 
T' G contains all the directed edges (ef,, e l w ). Analogously, let bx, . . . , bi be the black vertices incident to w in 
G and different from b, and let /* = {w, bi}. Then T' G contains all the directed edges {e w ,fl) (see Figure[8|). 
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Figure 8. Directed edges incident to et and e' w , e — {b, wi}, e' = {6, W2} 

Let us get back now to the expression of the Feynman function F(A(s)(x)) as the square root det(7 — 
A'(s)(x)). For example, one can apply the following formula (called Fermionic expansion): 

Lemma 2.21. Let A : V — » V be a linear operator acting on finite dimensional (complex) vector space 
V , dim(V) = size(A) = n. Let A' 1 A : A l V — > A' l V be the induced operator acting on A l V , the i-th skew 
symmetric power of V . Then 

n 

(19) 6eb(I-A) = ^2(-l) i Tr(A i A). 

i=0 

Proof. The functions det, Tr do not depend on the representative in the conjugacy class {A — > gAg~ x \g S 
Gl(V)}, so it is sufficient to consider the basis of V in which A has Jordan block structure. We shall stick to 
the generic case of semi-simple conjugacy class, i.e. all Jordan blocks have size 1. The proof in the presence 
of a non-trivial Jordan block follows from the generic case by continuity. 
Let Spec(A) = {\j}j=x,..., n - The left hand side of ((19)) reduces to 

n n 

det(I-A) = =£(-l) i «<(Ai > ... 1 A n ), 

i=0 i=0 

where Si (Aj., . . . , A„) is the homogeneity i symmetric function in variables Ai, . . . , A„. 

The right hand side of (fT§|) is evaluated as follows. Let {ei,...,e„} be the basis of V consisting of 
eigenvectors of A subordinate to the eigenvalues Ai, . . . , A„. The matrix A k A acts on the basis element 
e ix A ■ ■ ■ A e ik , 1 < i\ < i 2 < ■ ■ ■ < ik < n, by 

AM : e h A • • • A e ih -> Ae h A • • • A Ae lk = (A n . . . \ lk )e h A ■ ■ ■ A e ik . 
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The trace of /\ k A, computed in the basis {e^ A • • • A e.i k }i<i 1 <i 2 <—<ik<n 

of A k V, results in Tr(A k A) = 

Sfe(Ai, . . . , A n ) (note that dim{/\ k V) = (?)) and the alternating sum of Tr(/\ k A) over k G {0, 1, . . . , n} is 
equal to the left hand side. The proof follows. 

□ 

Let us mention a more general motivating geometrical and homological context for the Fermionic formula, 
see |rsl) . |rs2| . Let (M,g) be a compact Riemannian manifold and E — > M a vector bundle on M. We can 
associate to the spectrum of Laplace operator acting on i-forms twisted by E, 

A g : A l T*M ®E^ A i T*M ® E, 

its spectral zeta function 

A je Spec(A*)\{0} 

defined by analytic continuation of the right hand side convergent for Re(s) 3> to the whole complex plane 
C. Here Spec(A g ) = {Xj}j- The zeta function regularized determinant of is defined by 

(20) det'(Al) :=exp-^a s }l-«. 
Then the analytic torsion T(E, g, A g ) is 

dim(M) 

(21) T(E,g,Al):=( J[ {def A^^ 1 ^ . 

3=0 

The logarithm of analytic torsion is then half of the weighted sum of functional determinants for Laplace 
operators. In particular, observe that in the case of analytic torsion on a Riemann surface (X, g) we have 

(22) T a { 9l A^)^{det'A^)K 

where A®' 1 is the Laplace operator acting on (0, l)-forms of (X, g), rs2 . 

The discrete analogue of 1-forms was considered by Mercat in [m], see also Lovasz [lo]. A 1-form is a 
function on the set of directed edges of an embedded graph. This parallels to our matrix A'(s)(x), which is 
being indexed by the directed edges of the critically embedded graph G. 

It is clear that the previous zeta function regularization boils down to the determinant computation in 
the case of an operator T acting on a finite dimensional vector space. Let {Ai, . . . , A^} C Spec(T) be the 
subspace of non-zero eigenvalues of T. Its spectral zeta function is 

k 

Ct( S ) = Ev s 

and so 

^£>p)| s =o = -£ln(A,) = -HU.\i) = -Hdet'{T)). 

z=l 2=1 i—i 

In conclusion, exp - 3H T ( s )l s=0 is equal to the regularized determinant det'{T) of T. 

3. Appendix: Pfaffian method, Dimers at criticality and Determinants 

This paper studies a question whether one can define, using the geometric data provided by a critical 
embedding of a graph G in a Riemann surface, weights of the edges of G in a way that the matrix of these 
weights captures the basic properties of discrete Dirac operator and the partition function of the Dimer or 
Ising models with these weights may be explicitly evaluated. 

For the sake of completeness we briefly review the Pfaffian method which gives a general way to enumerate 
these partition functions for embedded graphs. 
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3.1. Kasteleyn orientations. The first step in this theory is the reduction of £(G,x) to V(Gg,x), where 
G$ is obtained from G by local changes which do not affect genus of the Riemann surface X in which G is 
embedded. This construction (see e.g. |lmj ) is not relevant here, so we omit it and concentrate on the Dimer 
partition function V{G,x). 

Assume the vertices of G are numbered from 1 to n. An orientation of G is obtained by prescribing one 
of the two possible directions to each edge of G. If D is an orientation of G, we denote by A(G, D) the 
skew-symmetric adjacency matrix of D defined as follows: The diagonal entries of A(G, D) are zero, and the 
off-diagonal entries are 

A(G, D)ij = ±x e , 

where the sum is taken over all edges e connecting vertices i and j, and the sign in front of x e is 1 if e is 
oriented from i to j in the orientation D, and —1 otherwise. It is well known that the Pfafhan of A(G, D) 
counts perfect matchings V(G) of the graph G with signs: 

Pfaf A(G,D) = sign(M,£>) J| x e , 

MeV{G) eEM 

where sign(M , D) = ±1. We use this as the definition of the sign of a perfect matching M with respect to 
the orientation D. The following statement is the basic result in the field: 

Theorem 7 (Kasteleyn [ka], Galluccio-Loebl [gl], Tesler [t], Cimasoni-Reshetikhin |SE])- If G embeds into a 
Riemann surface of genus g, then there exist 4 9 orientations Di (i = 1,...,4 9 ) of G such that the per- 
fect matching polynomial V(G, x) can be expressed as a linear combination of the Pfaffian polynomials 
Pfaf A{G, Di)(x). Moreover, the orientations Di are Kasteleyn orientations, i.e., they satisfy the following 
property: if F is a facial cycle of G then F has an odd number of edges oriented in Di in agreement with 
the clockwise traversal of F. 

The formula of Theorem [7] is called the Arf-invariant formula, as it is based on the property of the 
Arf invariant for quadratic forms in characteristic two. As far as we know, the relationship with the Arf 
invariant was first observed in [nrj . In |kuj , Kuperberg introduced a generalization of the notion of Kasteleyn 
orientations, called Kasteleyn flatness. 

3.2. Kasteleyn curvature. We recall that a graph G = (V, E) is called bipartite if the set of the vertices 
V may be partitioned into two sets W, B so |e H W\ = |e CI £?| = 1 for each e £ E. In this and the next one 
subsection we restrict ourselves to finite bipartite graphs G = (W,B,E), where W, B are the two edge-less 
sets of vertices of G and V — W U B. We call the vertices in W white and the vertices in B black and 
assume that G has at least one perfect matching. In particular, we restrict to the case of equal cardinalities 
\W\ = \B\. By a cycle in G we mean a subset of edges C C E(G) which form a cycle. The cycle C can be 
decorated by one of the two possible orientations, i.e. one of the two possible ways of going around C. By 
an oriented cycle we mean a cycle decorated with an orientation. In this subsection we also use O to denote 
the orientation of G in which each edge is oriented from its black vertex to its white vertex. 

Definition 3.1. Let G be given with the weights w(e),e € E(G). Let C be an oriented cycle of G. We 
define 

c(C) = ( _ 1) |c|/ 2+i n eg c + ^(e) 

as the Kasteleyn curvature of C. Here C + denotes the subset of edges of C whose orientation inherited from 
C coincides with their orientation in O, and C_ = C \ C + . 

Definition 3.2. Let G be a weighted graph, u>(e),e G E(G), embedded in X. We say that G is Kasteleyn 
flat if c(F) = 1 for each face F of the embedding, arbitrarily oriented. 

Let G = (W,B,E) be a bipartite graph equipped with the weights w(e),e e E(G). Let us fix a linear 
ordering on the set BL)W such that the elements of B precede the elements of W. This allows to introduce 
D(w) as (BUW) x (BUW) skew-symmetric matrix defined by D(w) uv = w{uv)x e if u e B and e = uv e E{G) 
resp. D(w) uv = — w(uv)x e if u £ W and e = uv 6 E(G), D(w) uv — otherwise. We further denote by 
Db{w) the (B x W)— block of D(w). For Al a perfect matching of G we denote by t(M) the coefficient of 
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IleeAf Xe m det(-Ds(tu)). Note that t(M) = sign(M) JleeM w ( e )> where sign(M) is the sign of M relative to 
the fixed linear ordering on the set of vertices B U W . 

The significance of flatness of the Kasteleyn curvature rests on the following observation (see |ku] ) . 

Proposition 3.3. Let C be a cycle of G embedded in X so that C is a symmetric difference of two perfect 
matchings M,N and the orientation of C coincides with the orientation of the edges of M in O. Then 

t(M) 



t(N) 

This result leads to 



c(C). 



Proposition 3.4. // G is embedded in the plane and its weight-function is Kasteleyn flat then t{M) is 
constant. 

We are able to prove a stronger statement. Let G be embedded in a Riemann surface X of genus g. We 
say that the weigh-function w is simple Kasteleyn flat if it is Kasteleyn flat and moreover w(e) G {1, —1} 
for each e G E. The function w may be viewed as assigning orientation O(w) obtained from orientation O 
by reversing the orientation of all edges with negative weight. If w is simple Kasteleyn flat the orientation 
D{w) is a Kasteleyn orientation, i.e. each face has an odd number of edges oriented clockwise. The next 
definition is motivated by [crj . 

Definition 3.5. We say that two weight-functions w,w' are equivalent if w can be obtained from w 1 by 
finite number of vertex multiplications, where a vertex multiplication consists in choosing a vertex v and a 
complex number c ^ together with multiplication by c the weights of all edges incident with v. 

Proposition 3.6. Let G be a graph embedded in a closed Riemann surface X of genus g and let w be a 
Kasteleyn fiat weight function, satisfying in addition c(C) G {1,-1} for each cycle C. Then w is equivalent 
to a simple Kasteleyn fiat weight-function w' , i.e., to a Kasteleyn orientation. 

Proof. We may assume that G is connected. Let T = (V(G), £"), E 1 C E(G) be a spanning tree of G. We 
can clearly perform vertex multiplications (in G) in a way the resulting weight function w' satisfies w'(e) = 1 
for each e € E'. Let e G E\E' . Then necessarily w'{e) G {1,-1} since e forms a cycle (say, denoted by C) 
with a subset of E' and c(C) has values in {1, —1}. □ 

Proposition 13.61 along with Theorem [7] immediately imply 

Corollary 3.7. Let G be a graph embedded in a closed Riemann surface X of genus g and w be a Kasteleyn 
flat weight function satisfying c(C) G {1, —1} for each cycle C. Then V{G,x) is a linear combination of A 9 
determinants of matrices obtained from Db{w) by multiplying some entries by — 1. 

3.3. Kasteleyn curvature and criticality. The computation of the Dimer partition function V(G,x) for 
a critically embedded bipartite graph G is based on the following strategy (see [k], jej): define the weight 
function w using discrete geometric information contained in the data of the critical embedding and use this 
information to prove that the kernel of the corresponding matrix D{w) has properties desired for a discrete 
analogue of the Dirac operator. Next, prove that w is Kasteleyn flat and insert it into Corollary 13. 71 

This approach works for planar bipartite graphs as well ( [k] ) , since in this case each Kasteleyn flat weight 
function satisfies c(C) G {1,-1} for each cycle C, see ku . On the other hand the assumptions of Corollary 
13.71 are quite restrictive for non-planar surfaces. By Proposition 13. 61 they are equivalent to w(e) G {1, —1} 
for each edge e G E. Moreover, it is shown in [nj that the theory of Kasteleyn flatness can not go beyond 
Corollary [XT] 
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